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ABSTRACT  INTRODUCTION 


The  wake  behind  a  two-bladed  model  rotor  in 
light  climb  was  measured  using  particle  image  ve- 
locimetry,  with  particular  emphasis  on  the  develop¬ 
ment  of  the  trailing  vortex  during  the  first  revolution 
of  the  rotor.  The  distribution  of  vorticity  was  dis¬ 
tinguished  from  the  slightly  elliptical  swirl  pattern. 
Peculiar  dynamics  within  the  “void”  region  may  ex¬ 
plain  why  the  peak  vorticity  appeared  to  shift  away 
from  the  center  as  the  vortex  aged,  suggesting  the  on¬ 
set  of  instability.  The  swirl  and  axial  velocities  (which 
reached  44  and  12  percent  of  the  rotor-tip  speed,  re¬ 
spectively)  were  found  to  be  asymmetric  relative  to 
the  vortex  center.  In  particular,  the  axial  flow  was 
composed  of  two  concentrated  zones  moving  in  oppo¬ 
site  directions.  The  radial  distribution  of  the  circu¬ 
lation  rapidly  increased  in  magnitude  until  reaching 
a  point  just  beyond  the  core  radius,  after  which  the 
rate  of  growth  decreased  significantly.  The  core-radius 
circulation  increased  slightly  with  wake  age,  but  the 
large-radius  circulation  appeared  to  remain  relatively 
constant.  The  radial  distributions  of  swirl  velocity 
and  vorticity  exhibit  self-similar  behaviors,  especially 
within  the  core.  The  diameter  of  the  vortex  core  was 
initially  about  10  percent  of  the  rotor-blade  chord, 
but  more  than  doubled  its  size  after  one  revolution  of 
the  rotor.  According  to  vortex  models  that  approx¬ 
imate  the  measured  data,  the  core-radius  circulation 
was  about  79  percent  of  the  large-radius  circulation, 
and  the  large-radius  circulation  was  about  67  percent 
of  the  maximum  bound  circulation  on  the  rotor  blade. 
On  average,  about  53  percent  of  the  maximum  bound 
circulation  resides  within  the  vortex  core  during  the 
first  revolution  of  the  rotor. 


Unlike  the  wake  of  a  fixed-wing  aircraft,  the  wake 
of  a  helicopter  rotor  follows  a  trajectory  that  repeat¬ 
edly  passes  near  portions  of  the  rotor  system  before 
it  is  finally  convected  downstream.  As  a  consequence, 
the  rotor  wake  not  only  has  a  considerable  influence 
on  the  performance  of  the  rotor,  but  it  is  principally 
responsible  for  the  high  levels  of  noise  and  vibration 
that  result  from  encounters  between  the  rotor  blades 
and  their  collective  wake  (ref.  1). 

Because  of  the  complexity  of  the  rotor  wake, 
especially  concerning  the  intensity  and  size  of  the 
dominant  trailing  vortex  and  its  increasingly  irreg¬ 
ular  trajectory,  computational  modeling  of  the  flow 
field  continues  to  be  a  challenge.  An  impressive  ef¬ 
fort  was  made  recently  to  calculate  the  performance 
and  aeroacoustics  of  the  rotor  without  knowing  the 
precise  structure  of  the  trailing  vortex  (ref.  2).  An 
equally  notable,  although  computationally  demand¬ 
ing,  computational  fluid  dynamics  (CFD)  code  has 
been  developed  that  attempts  to  accurately  capture 
the  details  of  the  rotor  wake  development  (ref.  3). 
To  achieve  this  level  of  accuracy,  numerical  dissipa¬ 
tion  must  be  minimized  by  imposing  very  small  spa¬ 
tial  and  temporal  step  sizes.  Even  though  this  highly 
resolved  computational  grid  is  selectively  structured 
in  the  flow  field  only  where  it  is  most  needed,  an  ex¬ 
traordinary  length  of  time  is  required  to  perform  the 
calculations.  However,  because  of  a  more  relaxed  grid 
requirement,  shorter  computational  times  have  been 
reported  using  a  CFD  formulation  that  incorporates 
a  vorticity-confinement  concept  (ref.  4).  Reliable  ex¬ 
perimental  measurements  are  needed  to  support  such 
efforts. 

As  new  experimental  techniques  emerge,  the  ro¬ 
tor  wake  is  frequently  revisited  in  order  to  extract 
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more  detailed  information  about  the  structure  of  the 
trailing  vortex  (refs.  5  -  18).  For  example,  not  only 
does  the  stereo  particle  image  velocimetry  (PIV)  mea¬ 
surement  technique  render  all  three  components  of  ve¬ 
locity  over  a  suflSciently  large  cross  section  of  the  flow, 
but  the  vector  fields  from  a  sequence  of  images  taken 
during  successive  revolutions  of  the  rotor  can  also  be 
conditionally  ensemble  averaged  to  remove  the  effects 
of  vortex  wander  (ref.  17).  Given  this  opportunity 
to  acquire  data  more  efficiently,  highly  resolved  (in 
space  and  time)  data  on  the  development  and  aging  of 
the  trailing  vortex  is  now  possible.  In  addition  to  en¬ 
hancing  our  physical  understanding  of  the  rotor  wake, 
specific  characteristics  (such  as  composition,  size,  in¬ 
tensity,  and  location)  of  the  trailing  vortex  at  different 
wake  ages  will  enable  CFD  models  to  be  validated  on 
a  more  detailed  level. 

Many  questions  have  been  raised  concerning  un¬ 
expected  features  appearing  within  the  trailing  vor¬ 
tex  (ref.  6),  the  composition  of  the  vortex  relative 
to  blade  loading  (ref.  2),  vortex  definitions  (ref.  19), 
opinions  dealing  with  the  development  and  demise  of 
vortex  systems  (ref.  20),  and  the  appropriateness  of 
simple  vortex  models  (ref.  21).  In  response  to  these 
questions,  an  experiment  was  performed  to  provide 
data  that  focuses  specifically  on  cross  sections  of  the 
trailing  vortex  during  a  complete  revolution  of  the  ro¬ 
tor.  The  data  are  of  sufficient  quality  to  permit  vari¬ 
ous  modeling  uncertainties  to  be  resolved.  This  paper 
describes  the  test  facility,  model,  flow  measurement 
technique,  and  the  data-acquisition  and  -processing 
procedure.  A  discussion  of  the  results  is  provided 
along  with  specific  findings. 


NOMENCLATURE 


a 

speed  of  sound,  meters  per  second  (m/s) 

A 

area  enclosed  by  circulation  path,  m^ 

c 

chord  length,  m 

Cq 

torque  coefficient,  torque/p7rf2^R® 

Ct 

thrust  coefficient,  thinst / pnQ^R^ 

dh 

distance  between  velocity  peaks,  horizontal  cut 

dy 

distance  between  velocity  peaks,  vertical  cut 

M 

hover-tip  Mach  number,  RQ./ a 

N 

number  of  points  to  be  averaged 

Nt 

number  of  rotor  blades 

R 

radius  of  rotor-tip  path,  m 

Re 

hover-tip  Reynolds  number,  cRQ,/v 

r,0,z 

coordinates  relative  to  rotor  axis 

u,t;,  w 

velocity  components,  m/s 

UcyWc 

vortex-convection  velocities,  m/s 

Ur^Wr 

velocities  relative  to  UcyWc,  m/s 

a 

angle  of  position  vector,  degrees  (deg) 

p 

angle  of  u  -f-  u?  vector,  deg 

r 

circulation,  m^/s 

To 

circulation  at  core  radius,  m^/s 

r. 

total  vortex  circulation  in  far  field,  m^/s 

5 

turbulent /laminar  viscosity,  e/z^ 

€ 

eddy  viscosity,  m^/s 

c 

radius  measured  from  void  center,  mm 

Co 

radius  of  vortex  core,  millimeters  (mm) 

CcO 

value  of  Cc  at  trailing  edge,  mm 

c 

C/Cc 

0 

local  angle  of  rotor  blade,  deg 

K 

radius  measured  from  vorticity  centroid,  mm 

X 

dimensionless  hover-inflow  velocity,  Wi/Q>R 

u 

kinematic  molecular  viscosity,  m^/s 

P 

fluid  density,  kg/m^ 

a 

rotor  solidity,  Niyc/i^R 

9 

swirl  velocity,  m/s 

9c 

swirl  velocity  at  core  radius,  m/s 

9 

9l9c 

azimuthal  wake  age  measured  from  blade,  deg 

vorticity  normal  to  r  —  z  plane,  per  second 

Q 

rotor  revolutions  per  minute,  rpm 

ABBREVIATIONS 

CCD 

Charged  coupled  device 

CFD 

Computational  fluid  dynamics 

LDV 

Laser  Doppler  velocimetry 

PIV 

Particle  image  velocimetry 

TEST  DESCRIPTION 
Test  Facility 

This  experiment  was  performed  in  the  settling 
chamber  of  the  7-  by  10-Foot  Subsonic  Wind  Tunnel 
at  Ames  Research  Center  (fig.  1).  This  location  was 
selected  to  take  advantage  of  the  relatively  large  cross 
section  of  the  tunnel  (9.4  m  x  10.0  m)  and  to  better 
isolate  the  flow  field  around  the  rotor  from  ambient 
wind  fluctuations  external  to  the  facility. 

A  model  rotor  was  placed  in  the  center  of  the  set¬ 
tling  chamber  and  mounted  horizontally  with  its  axis 
of  rotation  coincident  with  the  centerline  of  the  test 
section.  Aside  from  the  influence  of  the  tunnel  walls 
(approximately  4.6  m  away)  and  the  contraction  into 
the  test  section  (with  an  area  ratio  of  14:1),  the  ro¬ 
tor  wake  was  allowed  to  pass  through  the  test  section 
without  interference.  A  curtain  was  drawn  across  the 
air-exchanger  section  to  prevent  the  wake  from  re¬ 
circulating  around  the  tunnel  and  contributing  to  an 
unwanted  climb  condition.  In  this  configuration,  any 
portion  of  the  wake  reaching  this  section  would  be 
expelled  through  the  air-exchanger  doors  and  fresh 
stagnant  air  would  be  drawn  by  the  rotor  into  the 
tunnel  circuit  (as  required  by  mass  continuity). 
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Video  cameras  were  mounted  upstream  and  to 
the  side  of  the  rotor  to  monitor  critical  elements  dur¬ 
ing  the  test  (fig.  2),  It  is  customary  to  infer  the 
velocity  in  the  test  section  from  the  differential  pres¬ 
sure  between  the  settling  chamber  and  the  test  sec¬ 
tion.  However,  this  procedure  does  not  reveal  the 
nonuniform  character  of  the  flow  across  the  test  sec¬ 
tion  (due  to  blockage  by  the  PIV  apparatus  in  the 
settling  chamber)  and  is  not  suflSciently  accurate  at 
low  speeds.  Therefore,  the  velocity  in  the  test  sec¬ 
tion  due  to  the  rotor  wake  was  measured  with  a  vane 
anemometer  (1  to  39.6  m/s  range)  attached  to  an  ex¬ 
isting  tunnel  traverse.  In  the  settling  chamber,  where 
the  velocity  was  considerably  lower,  a  hot-film  probe 
was  used  to  measure  velocity  (0  to  2.4  m/s  range)  di¬ 
rectly  upstream  of  the  rotor  at  a  distance  of  about  8.5 
meters  from  the  tip-path  plane. 

Model  Rotor 

The  rotor  used  in  this  experiment  was  a  1/7-scale 
model  of  the  full-scale  AH-l/OLS  rotor.  The  profile 
of  this  blade  was  based  on  a  simple  trailing-edge  ex¬ 
tension  to  the  shape  of  a  NACA  0012  airfoil  (fig.  3). 
The  rotor  consisted  of  two  blades,  each  with  a  chord 
of  103.9  mm,  maximum  thickness  of  10.2  mm,  and  a 
span  of  957.6  mm  measured  from  the  axis  of  rotation 
to  the  blade  tip  (fig.  4).  Beginning  at  25-percent  ra¬ 
dius,  the  twist  of  each  blade  was  — 0.106^/centimeter 
(cm),  so  that  the  total  variation  in  geometric  twist 
along  the  span  was  7.55^  (from  the  center  of  rota¬ 
tion  to  the  rotor  tip).  Collective  pitch  settings  were 
referenced  to  the  geometric  pitch  angle  at  75-percent 
radius.  Therefore,  a  collective  pitch  setting  of  zero 
meant  that  the  blade  pitch  was  about  5.0®  at  the  root 
(held  fixed  for  locations  less  than  R/A)  and  —2.5®  at 
the  tip.  The  rotor  hub  had  a  built-in  blade  precone 
of  1.30®. 

Information  related  to  safety  of  flight  was  ob¬ 
tained  by  flap-  and  chord-bending  gages  mounted  on 
the  blade.  The  drag  link  used  to  adjust  the  blade 
sweep  was  instrumented  with  a  strain  gage  that  re¬ 
sponded  to  the  aerodynamic  drag  on  the  rotor,  offer¬ 
ing  a  redundant  measurement  of  chord- wise  bending. 
Depending  on  the  sweep  of  the  blade  and  the  loca¬ 
tion  of  the  center  of  mass  of  the  blade,  this  gage  could 
also  be  loaded  by  a  centrifugal  force  induced  moment. 
Strain  gages  were  placed  on  the  pitch  horns  that  trans¬ 
mitted  commanded  pitch  to  the  rotor  blades.  Reac¬ 
tions  to  these  commands  resulted  in  bending  loads  on 
these  elements  that  were  proportional  to  the  torsional 
moments  on  the  blades.  A  transducer  for  measuring 
the  angular  tilt  of  the  rotor  (called  the  teeter  pot)  was 
mounted  on  top  of  the  hub.  The  rotor  operator  mon¬ 
itored  the  signals  from  each  of  these  gages.  Because 
the  center  of  mass  of  each  rotor  blade  was  located  aft 


of  the  feathering  axis  (aligned  with  the  quarter  chord), 
it  was  necessary  to  sweep  the  blades  forward  1.20®  to 
minimize  steady  chordwise  bending  due  to  centripetal 
force.  The  blades  were  strobed  at  2/revolution  (rev) 
to  verify  that  the  blades  were  tracking.  A  blade  sweep 
of  1.20®  corresponded  to  a  0.19c  forward  movement  of 
the  tip. 

The  model  rotor  was  attached  to  a  teetering  hub, 
which  in  turn  was  mounted  on  the  VSB-54  force  and 
moment  balance.  The  most  significant  element  of 
the  balance  in  the  present  test  was  the  axial  gage, 
which  had  a  capacity  of  159  kilograms  (kg)  and  pro¬ 
vided  a  measure  of  the  rotor  thrust.  Rotor  torque 
was  measured  using  torsion  gages  mounted  on  a  flex 
coupling  located  between  the  balance  and  the  trans¬ 
mission,  and  had  a  capacity  of  74.7  kg-m.  The  rotor 
was  driven  by  a  90-horsepower  (HP)  electric  motor 
through  a  4:1  transmission.  The  complete  rotor  and 
drive  system  (fig.  5)  was  supported  on  a  tripod-like 
structure  that  wets  rigidly  attached  to  the  frame  of  the 
wind  tunnel.  This  setup  placed  the  tip-path  plane  of 
the  rotor  about  1.8  m  upstream  of  the  beginning  of 
the  contraction  section  of  the  tunnel,  well  within  the 
protective  steel- wall  section  of  the  tunnel. 

Stereo  PIV  Setup 

The  primary  elements  of  the  PIV  system  con¬ 
sisted  of  a  pulsed  laser,  optics  to  spread  the  beam  into 
a  thin  sheet,  charged  coupled  device  (CCD)  cameras 
to  record  the  images,  and  a  calibration  target.  The 
plane  of  the  light  sheet  was  horizontal  and  passed  just 
below  the  trailing  edge  of  the  rotor  when  the  rotor  was 
at  zero  azimuth  (fig.  6).  In  order  for  both  cameras  to 
receive  close  to  the  same  amount  of  forward-scattered 
light,  two  opposing,  coplanar  light  sheets  were  estab¬ 
lished  (figs.  6  and  7).  To  avoid  interference  with  the 
model  supports,  the  cameras  were  placed  at  different 
distances  away  from  the  light  sheet  (nominally  1.5  m) 
and  at  different  angles.  Two  Nd:YAG  lasers,  120  mil- 
lijoules  (mJ)  each,  were  positioned  and  synchronized. 
Cameras  and  lasers  were  triggered  off  of  a  1/rev  ref¬ 
erence  signal  that  was  keyed  to  the  rotor.  This  sig¬ 
nal  was  obtained  by  directing  a  small  continuous- wave 
laser  beam  onto  a  small  mirror  attached  to  the  hub 
of  the  rotor  and  intercepting  the  reflected  beam  with 
a  fast-response  photodetector.  Acquisition  software 
was  used  to  delay  the  reference  signal  in  order  to  ob¬ 
tain  data  at  different  wake  ages. 

Images  were  acquired  with  8-bit  CCD  cross¬ 
correlation  cameras,  each  having  a  sensor  array  of  2K 
X  2K  pixels  (nominal).  Each  pixel  measured  7.4  mi¬ 
crons  (/im)  on  a  side.  The  cameras  were  operated 
in  a  double- exposure  mode  to  acquire  two  noninter¬ 
laced,  full-frame  images  during  a  single  frame  interval. 
These  were  nonstandard  video  cameras  that  could  be 
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externally  triggered  and  driven  at  any  frequency  up  to 
15  Hertz  (Hz)  in  double-exposure  mode.  A  computer 
interface  provided  control  over  gain,  contrast,  black 
level,  and  trigger  mode.  Because  the  cameras  were 
located  at  different  distances  from  the  object  plane, 
one  used  a  55-mmlens,  while  the  other  used  an  85-mm 
lens.  The  lenses  were  remotely  translated  to  focus  on 
the  centerline  of  the  image  area,  and  the  sensor  (lo¬ 
cated  inside  the  camera  body)  was  rotated  about  its 
centerline  to  satisfy  the  Scheimpflug  condition^. 

Each  NdrYAG  laser  was  frequency  doubled  to 
provide  a  beam  having  a  wavelength  of  532  nanome¬ 
ters  (nm),  a  pulse  width  of  about  4  nanoseconds  (ns), 
and  power  of  120  mJ.  The  laser  beam  initially  had  a 
diameter  of  about  5  mm  and  a  divergence  of  0.60  mil- 
liradians  (mrad).  In  this  test,  each  laser-light  sheet 
was  formed  by  first  passing  the  beam  through  a  set 
of  spherical  lenses  (consisting  of  focal  lengths  of  —500 
mm  and  +500  mm).  By  adjusting  the  distance  be¬ 
tween  these  two  lenses,  the  beam  waist  (or  sheet  thick¬ 
ness,  approximately  1  mm  thick)  could  be  positioned 
in  the  region  imaged  by  the  cameras.  Two  cylindri¬ 
cal  lenses  (with  focal  lengths  of  —12.7  mm  and  +25.4 
mm)  were  used  to  expand  the  beam  in  the  horizon¬ 
tal  direction.  The  distance  between  the  sheet-forming 
optics  and  the  imaged  area  was  about  1.5  m. 

The  lasers,  sheet-forming  optics,  and  cameras 
were  rigidly  mounted  on  a  common  structure,  which 
could  be  traversed  in  three  orthogonal  directions. 
This  arrangement  allowed  the  cameras  to  be  focused 
only  on  the  trailing  vortex  region  of  the  flow,  and  then 
traversed  as  required  to  follow  the  position  of  the  vor¬ 
tex  at  different  wake  ages. 

Prior  to  taking  data,  a  flat  aluminum  plate,  with 
a  15-  by  23-cm  rectangle  precisely  milled  into  the  sur¬ 
face,  was  used  to  calibrate  the  camera  views.  The 
calibration  plate  was  leveled  and  the  light  sheets  ad¬ 
justed  so  as  to  evenly  graze  the  surface.  The  lens  of 
each  camera  was  traversed  along  its  optical  axis  and 
the  CCD  sensors  were  rotated  until  a  uniform  focus  on 
the  calibration  target  was  achieved.  The  calibration 
supplied  all  the  parameters  required  to  accurately  de- 
warp  the  images  and  construct  the  three-component 
velocity  fields. 

Proper  seeding  of  the  flow  is  critical  for  obtaining 
accurate  PIV  measurements.  The  seed  particles  must 
be  evenly  distributed  and  of  sufficient  density  to  de¬ 
fine  the  flow  without  altering  its  physical  properties. 
The  particles  must  also  be  small  enough  to  accurately 
follow  the  flow  (especially  challenging  in  accelerating 
flows),  yet  large  enough  so  that  they  scatter  a  suf¬ 
ficient  amount  of  light  to  be  detected.  The  particle 

^Uniform  focus  requires  that  the  light  sheet,  camera  lens, 
and  camera  sensor  planes  intersect  along  a  common  line  (see 
fig.  7). 


generator  used  in  this  test  employed  an  inert  gas  to 
atomize  a  nontoxic,  pharmaceutical-grade  mineral  oil. 
The  mist  was  vaporized  and  then  condensed  before 
being  released  into  the  flow.  The  particle  size  was  es¬ 
timated  by  the  manufacturer  to  be  less  than  0.5  pm. 
The  particle  generator  was  located  upstream  of  the 
settling  chamber  so  that  the  particles  would  mix  with 
the  air  entering  the  tunnel  circuit  before  being  drawn 
into  the  wake  of  the  rotor. 

Calibration,  data  acquisition,  and  image  process¬ 
ing  were  all  performed  using  Provision  software  devel¬ 
oped  by  Integrated  Design  Tools  (IDT).  The  princi¬ 
ples  of  the  PIV  technique  and  the  application  of  this 
software  to  the  study  of  vortex-dominated  flows  are 
discussed  in  greater  detail  in  reference  17. 

PROCESSING  WAKE  IMAGES 

After  completion  of  the  test,  the  data  were  pro¬ 
cessed  on  a  PC  workstation.  To  obtain  good  resolu¬ 
tion  of  the  trailing  vortex,  a  calculation  grid  measur¬ 
ing  92  horizontal  nodes  by  79  vertical  nodes  was  con¬ 
structed  over  the  region  of  significant  interest  in  the 
flow  field.  The  area  covered  by  the  grid  resulted  in 
approximately  12  pixels  between  nodes  in  both  direc¬ 
tions,  as  determined  by  (Grid  range  in  terms  of  pix- 
els)/(Number  of  nodes  —1).  The  interrogation  win¬ 
dow  was  set  at  22  pixels  on  a  side,  giving  a  45-percent 
overlap  according  to  (Interrogation  window  —  Interval 
between  nodes) /(Interrogation  window).  The  physi¬ 
cal  area  measured  in  each  image  zone  was  168  mm 
(with  an  axial  increment  of  1.8  mm)  by  106  mm  (with 
a  radial  increment  of  1.4  mm).  The  1.8-  by  1.4-mm 
increments,  which  defined  the  calculation  step  sizes  in 
the  axial  and  radial  directions,  also  defined  the  reso¬ 
lution  of  the  measurements.  Therefore,  the  interroga¬ 
tion  window  (containing  484  pixels)  covered  an  area 
measuring  2.5  mm  in  the  radial  direction  by  3.4  mm 
in  the  axial  direction. 

Consider  an  interrogation  area  (fig.  8)  in  which 
10  particle  pairs  (the  minimum  number  required  to 
yield  a  valid  resultant  vector)  have  been  registered 
on  a  22-  by  22-pixel  image  of  a  small  region  in  a 
swirling  flow.  The  software  derives  a  representative 
in-plane  velocity  vector  based  on  a  second-order  cal¬ 
culation  of  all  the  individual  particle  displacements 
within  that  interrogation  area.  Note  that  the  effec¬ 
tive  “probe  size’^  is,  therefore,  less  than  the  interro¬ 
gation  area.  All  the  rotational  velocity  components 
within  the  vortex  are  accounted  for,  but  they  are  rep¬ 
resented  in  a  collective  sense  over  contiguous  regions 
defined  by  the  size  of  the  interrogation  area.  For  an 
initial  core  diameter  of  10  mm,  the  core  would  have 
been  covered  by  about  nine  interrogation  areas  (78.5 
mm^/8.5  mm^).  Considering  an  interrogation  overlap 
of  47  percent  means  that  about  18  vectors  would  be 
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calculated  within  the  vortex  core.  Given  that  the  size 
of  the  vortex  grows  as  it  ages,  the  number  of  vectors 
calculated  within  the  core  would  also  increase.  In  this 
present  test,  for  example,  a  vortex  that  had  grown  to 
a  diameter  of  20  mm  would  contain  about  72  vectors. 

A  single  file  for  each  wake  age  included  a  flag  for 
every  velocity  measurement  that  indicated  the  nature 
of  that  vector  (such  as  valid,  invalid,  interpolated,  re¬ 
calculated,  or  not  calculated).  Reduced  datafiles  con¬ 
taining  the  coordinates  of  each  calculation  node  and 
the  three  components  of  velocity  (in  terms  of  displace¬ 
ment)  were  stored  in  ASCII  format  and  then  trans¬ 
ferred  to  a  mainframe  computer  for  analysis.  The  ve¬ 
locity  components  were  first  converted  from  displace¬ 
ment  units  to  velocity  units  based  on  the  pulse  du¬ 
ration  for  that  particular  measurement.  The  physical 
coordinates  were  transformed  so  that  r  =  =  0  would 

correspond  to  the  center  of  the  rotor  hub  (fig.  6). 

Because  of  centripetal  effects,  heavier  seed  par¬ 
ticles  appear  to  have  been  driven  outward  from  the 
center  of  the  vortex  as  it  moved  through  the  fluid.  Be¬ 
cause  the  size  of  the  remaining  lighter  particles  were 
most  likely  much  smaller,  less  light  would  have  been 
scattered  from  this  region,  thereby  giving  it  the  ap¬ 
pearance  of  a  “void.”  Nevertheless,  the  particles  that 
were  present  inside  the  void  were  sufficient  to  obtain 
PIV  measurements.  In  searching  for  vortex  character¬ 
istics  that  can  be  known  with  certainty,  the  size  of  the 
“apparent  void”  cannot  always  be  equated  to  the  size 
of  the  core,  nor  can  the  center  of  the  void  be  assumed 
to  be  coincident  with  the  center  of  swirl.  On  the  other 
hand,  the  physical  location  of  the  void  was  an  indis¬ 
putable  feature  and  was,  therefore,  judged  to  be  the 
most  reliable  indicator  of  how  the  vortex  was  convect- 
ing  through  the  fluid.  In  addition,  in  order  to  account 
for  the  ever-present  vortex  meander,  the  coordinates 
of  the  voids  for  a  set  of  images  acquired  over  time  at 
a  particular  azimuth  behind  the  rotor  provided  a  very 
accurate  basis  for  conditionally  averaging  the  images. 

The  distribution  of  particles  around  an  apparent 
vortex  void  often  exhibited  a  peculiar  pattern  (fig.  9). 
Particles  were  not  always  uniformly  distributed  in  the 
flow  as  it  approached  the  rotor.  In  these  instances, 
some  form  of  particle  “banding”  frequently  occurred 
adjacent  to  the  void  region.  Calculations  have  been 
made  that  show  that  banding  is  also  a  natural  con¬ 
sequence  of  the  balance  of  forces  (drag,  centrifugal, 
and  Coriolis)  that  act  on  the  seed  particles  within  a 
vortex  (ref.  22).  Small  particles  settle  in  an  annular 
region  that  is  smaller  than  the  vortex  core,  compared 
to  larger  particles  that  yield  a  band  that  is  larger  than 
the  vortex  core.  Regardless  of  how  uniform  the  flow 
appeared  to  be  seeded  in  the  present  experiment,  a 
fairly  broad  range  of  void  sizes  was  observed  at  all 
wake  ages  (fig.  10).  The  average  size  of  the  void  in¬ 


creased  rapidly  during  the  first  100"^  of  azimuth,  after 
which  the  rate  of  increase  diminished  substantially. 

In  this  study,  the  sequence  for  extracting  the 
characteristics  of  the  trailing  vortex  began  with  a 
physical  examination  of  each  raw  image  to  deter¬ 
mine  the  size  and  location  (relative  to  a  superimposed 
mesh)  of  the  apparent  void  (fig.  11).  At  each  wake 
age  the  images  were  averaged  to  smooth  out  the  small 
irregularities  that  were  invariably  present  in  the  flow. 
Because  the  vortex  was  located  at  different  positions 
from  image  to  image  (vortex  meander),  performing 
a  simple  average  would  have  risked  smearing  out  im¬ 
portant  details  concerning  the  structure  of  the  vortex. 
Instead,  the  data  were  averaged  only  after  aligning  the 
set  of  images  based  on  the  most  dominant  recurring 
feature  in  the  flow,  which  in  this  case  was  the  appar¬ 
ent  void. 

This  procedure  is  referred  to  as  conditional  en¬ 
semble  averaging,  the  condition  here  being  the  align¬ 
ment  of  all  the  images  in  the  set  based  on  their  void 
centers.  The  image  with  a  void  center  closest  to  the 
mean  location  that  was  established  for  the  set  of  im¬ 
ages  was  selected  as  the  “anchor.”  The  indices  of  all 
image  matrices  (that  is,  grid  points  i,j)  for  each  wake 
age  were  then  adjusted  according  to  the  offset  of  each 
void  from  the  coordinates  of  the  void  in  the  anchor 
image.  All  data  in  the  anchor  image  were  retained. 
However,  some  portion  of  all  other  images  with  ad¬ 
justed  indices  that  fall  outside  the  boundaries  of  the 
anchored  image  were  necessarily  discarded.  Hence,  as 
a  consequence  of  vortex  wander,  the  population  con¬ 
tributing  to  the  conditional  average  was  greatest  over 
the  interior  of  the  matrix  (or  calculation  grid) . 

At  this  intermediate  point,  the  conditionally  av¬ 
eraged  flow  field  at  each  wake  age  embodies  the  trail¬ 
ing  vortex  within  a  highly  nonuniform  downwash. 
The  time-averaged  downwash  for  a  single-bladed  rotor 
can  be  determined  by  averaging  the  flow  field  over  an 
entire  revolution  of  the  rotor  (refs.  8,  10).  When  the 
time-averaged  downwash  velocity  field  is  subtracted 
from  the  phase- averaged  velocity  field,  the  resulting 
velocity  profiles  tend  to  appear  more  symmetric  along 
the  extremities  of  the  vortex,  but  not  near  the  core 
radius  (ref.  8).  In  the  present  experiment,  this  pro¬ 
cedure  could  not  be  implemented  because  data  were 
not  obtained  over  the  required  half  revolution  in  each 
zone.  However,  the  axial  and  radial  locations  of  the 
voids  at  each  wake  age,  along  with  the  known  time 
intervals  between  each  wake  age,  were  known  and  did 
yield  the  convection  velocity  of  the  vortex  at  each 
wake  age.  These  components  of  the  convection  ve¬ 
locity,  Uc  and  Wc,  were  subtracted  from  the  phase- 
averaged  flow  field.  However,  because  the  axial  com¬ 
ponent  of  the  downwash  over  the  inboard  region  was 
significantly  higher  than  the  axial  convection  velocity 
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of  the  vortex  (i^c))  velocity  profiles  across  the  center 
of  the  vortex  cannot  be  expected  to  be  symmetric. 

To  extract  information  about  the  geometry  of  the 
vortex  (after  performing  averages  based  on  the  void 
centers  at  each  wake  age),  the  center  of  the  vortex 
was  assumed  to  be  at  the  center  of  swirl^.  Using  the 
void  center  as  a  starting  point,  the  surrounding  lo¬ 
cations  were  interrogated  to  determine  the  best  node 
for  which  the  sum  of  the  dot  products  of  two  unit 
vectors  produced  a  minimum  over  a  neighborhood^  of 
locations  surrounding  the  candidate  node.  The  coor¬ 
dinates  of  the  neighboring  node  relative  to  the  candi¬ 
date  node  defined  one  of  the  unit  vectors.  The  veloc¬ 
ity  vector  at  the  neighboring  node  defined  the  other 
unit  vector.  With  the  angles  of  these  two  unit  vectors 
denoted  by  a  and  (3  (fig.  12),  the  node  nearest  the 
center  of  swirl  was  said  to  be  found  when  the  following 
expression  was  a  minimum: 

N 

1]  |cos(a-/?)|/iV  (1) 

n=l 

where  N  is  the  total  number  of  neighboring  points 
considered.  Relative  to  the  center  of  swirl  for  a  pure 
vortex  in  a  stationary  flow,  a  and  /?  will  be  orthogonal 
for  all  neighboring  points  and  the  sum  will  be  exactly 
zero.  For  most  real  flows  the  sum  will  not  be  precisely 
zero,  and  the  procedure  works  best  when  the  neigh¬ 
borhood  under  consideration  does  not  extend  beyond 
about  one  or  two  core  diameters  (which  is  normally 
sufficient  to  cover  most  of  the  vorticity  in  a  single 
trailing  vortex). 

The  size  of  the  vortex  core  was  found  by  first  di¬ 
viding  the  surrounding  region  into  annular  zones  and 
calculating  the  average  swirl  velocity  and  the  associ¬ 
ated  mean  radius  of  each  annulus.  These  data  were 
fit  with  a  least-squares  spline  subject  to  the  condition 
that  the  resulting  curve  have  only  one  inflection  along 
its  inner  extent  (nominally  set  to  twice  the  estimated 
core  radius).  The  point  of  inflection  was  used  to  define 
the  radius  of  the  vortex  core,  because  that  is  where  a 
maximum  value  for  the  swirl  velocity  is  reached.  This 
procedure  was  considered  to  offer  a  more  rational  ap¬ 
proach  for  determining  the  size  of  the  vortex,  espe¬ 
cially  when  the  velocity  peaks  that  are  characteristic 
of  vortices  (and  upon  which  core  sizes  are  tradition¬ 
ally  based)  are  dependent  on  how  the  vortex  is  sliced 
(vertical,  horizontal,  or  otherwise). 

Vorticity  was  calculated  by  integration  (circula¬ 
tion  box  method)  according  to: 

^Neither  the  void  center  nor  the  swirl  center  wzis  necessarily 
at  the  center  of  the  vorticity  that  was  accreted  from  the  ro¬ 
tor  wake,  itself  having  the  form  of  a  spiraling  sheet  containing 
varying  amounts  of  vorticity. 

®  Because  the  quality  of  the  data  within  the  void  may  be  in 
question,  these  vectors  were  omitted  during  this  operation. 


ug  =  A  ^  ^  ^  (u,  iB)  •  d(r,  (2) 

where  the  direction  of  integration  was  such  that  the 
enclosed  area,  A,  was  on  the  left  of  the  integration 
path.  The  size  of  each  circulation  box  was  defined 
by  the  local  grid  spacing,  which  on  average  rendered 
A  =  2.5  mm2. 

Test  Conditions 

The  experiment  was  performed  at  a  constant  ro¬ 
tor  speed  of  870  rpm  (14.5  Hz),  which  for  a  rotor 
diameter  of  191.5  cm  corresponds  to  a  tip  speed  of 
87.2  m/s.  Based  on  an  average  ambient  temperature 
of  16.7  degrees  Centigrade  (®C),  the  Reynolds  num¬ 
ber  based  on  rotor-tip  chord  was  Re  =  0.6  x  10®  and 
the  Mach  number  was  M  =  0.26.  The  collective  pitch 
angle  of  the  rotor  was  fixed  at  IV. 

To  obtain  high-resolution  PIV  images,  the  cam¬ 
eras  were  focused  on  a  small  region  of  the  wake  behind 
the  trailing  edge  of  the  rotor.  The  camera  view  was 
radially  and  axially  translated  to  provide  three  zones 
of  coverage  (A,  E,  and  F).  The  PIV  measurement  lo¬ 
cations  (which  were  selected  to  follow  the  vortex  de¬ 
velopment  at  different  wake  ages),  along  with  the  laser 
pulse  delays  needed  to  capture  the  highest  velocities, 
are  shown  in  table  1.  Results  presented  later  in  this 
report  show  that  even  for  the  longest  pulse  delay,  5 
microseconds  (^s),  the  vortex  moved  less  than  1/1000 
of  the  core  radius  between  pulses.  The  blade  was  con¬ 
sidered  to  be  at  zero  azimuth  when  the  trailing  edge  of 
the  rotor  blade  was  horizontal.  When  the  blade  is  in 
this  position,  the  light  sheets  (both  of  which  intersect 
in  a  horizontal  plane)  pass  just  beneath  the  trailing 
edge  of  the  blade.  The  procedure  for  taking  images  at 
the  different  wake  ages  was  accomplished  by  delaying 
the  1/rev  signal  from  the  rotor,  thereby  allowing  the 
tip  of  the  rotor  blade  to  move  past  the  measurement 
plane  by  the  desired  azimuthal  increment.  The  data 
were  taken  at  a  nominal  increment  of  15^  (fig.  13), 
with  some  exceptions  as  required  to  avoid  reflections 
of  the  laser  light  from  any  portion  of  the  blade  surface 
that  might  damage  the  camera  sensors. 

RESULTS  AND  DISCUSSION 
Rotor  Performance 

Rotor  loads  and  tunnel  velocity  measurements 
(fig.  14)  were  taken  from  an  earlier  experiment  (ref. 
18).  The  colors  used  in  the  contour  plots  of  the  test- 
section  velocity  for  each  collective  angle  are  scaled  so 
that  red  represents  the  maximum  velocity  (and  blue 
the  minimum)  for  that  case.  Strictly  speaking,  the 
rotor  is  not  in  hover,  as  evidenced  by  the  magnitude 
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Table  1:  Measurement  locations 


Table  2:  Rotor  loads 


Wake  age,  deg 

Zone 

Pulse  delay,  ps 

0 

A 

5.0 

15 

A 

5.0 

30 

A 

5.0 

45 

A 

5.5 

60 

A 

5.0 

75 

A 

5.0 

90 

A 

4.5 

105 

A 

4.5 

120 

A 

4.5 

135 

A 

4.5 

150 

A 

4.5 

180 

E 

5.0 

195 

E 

3.5 

210 

E 

4.0 

225 

E 

4.0 

240 

E 

4.0 

255 

E 

4.0 

270 

F 

5.0 

285 

F 

5.0 

300 

F 

5.0 

315 

F 

5.0 

328 

F 

5.0 

344 

F 

5.0 

375 

F 

5.0 

390 

F 

5.0 

of  the  velocity  upstream  of  the  rotor.  Because  the 
flow  around  the  rotor  is  bounded  by  the  wind  tunnel 
walls,  and  because  the  rotor  wake  lowers  the  pressure 
in  the  test  section  (in  turn  inducing  a  finite  flow  ev¬ 
erywhere  ahead  of  the  rotor),  a  slight  climb  condition 
cannot  be  avoided  (see  ref.  18  for  a  full  discussion). 
Assuming  the  upstream  velocity  entirely  represents 
the  climb  velocity,  then  the  climb  rate  for  this  part  of 
the  experiment  can  be  calculated  from  Woo/R^,  and 
is  included  with  other  relevant  quantities  in  table  2. 

Rotor  Wake 

The  center  of  the  vortex  void,  the  center  of  vortic- 
ity,  the  center  of  swirl,  and  the  center  between  velocity 
peaks  are  all  considered  to  be  distinct  quantities.  For 
each  wake  age,  starting  with  ^  =  0^,  the  flow  field  has 
been  summarized  in  terms  of  quantities  that  focus  on 
the  trailing  vortex  (z.e.,  fig.  15).  Vortex  wander  was 
qualitatively  identified  by  plotting  the  locations  of  the 
vortex  voids  for  each  image  set  relative  to  the  calcu¬ 
lation  grid  that  was  used  to  compute  the  flow  field. 
Only  two  grid  steps  in  each  direction  are  shown.  To 
the  right  of  this  plot  is  a  sketch  of  the  rotor  show¬ 
ing  the  location  of  the  rectangular  zone  covered  by 


Rotational  speed  of  rotor,  rpm 

870 

Collective  pitch  angle,  deg 

11 

Solidity 

0.0691 

Torque,  foot-pounds  (ft-lbs) 

11.69 

Thrust,  lb 

39.76 

Upstream  velocity,  ft/s 

1.28 

Climb  rate 

0.0045 

Hover  figure  of  merit 

0.61 

Ct 

0.00657 

Cq 

0.00062 

Ct  1  O' 

0.095 

the  entire  calculation  grid  (92  steps  in  the  2:  direc¬ 
tion  and  79  steps  in  the  r  direction).  The  three  field 
plots  in  the  left-most  column  show  the  out-of-plane 
vorticity,  Qq  (contours  with  red  and  blue  hues  to  indi¬ 
cate  direction  and  a  continuous  range  of  saturations 
to  indicate  magnitude),  inplane  velocity  relative  to 
the  convection  velocity  of  the  vortex,  (vectors 

with  a  continuous  range  of  hues  to  indicate  magni¬ 
tude)  ,  and  the  out-of-plane  velocity,  v  (contours  with 
red  and  blue  hues  to  indicate  direction  and  a  continu¬ 
ous  range  of  saturations  to  indicate  magnitude).  The 
three  profiles  in  the  center  column  show  ,  and  v 

along  z  =  Constant  and  passing  through  the  center  of 
swirl.  The  central  portion  of  these  profiles  is  shown 
in  red  and  demarcated  with  dashed  vertical  lines  in 
order  to  identify  the  void  region.  The  horizontal  and 
vertical  dashed  lines  included  in  each  of  the  field  and 
profile  plots  mark  the  location  of  the  swirl  center. 

The  lower  plot  in  the  right-most  column  shows 
the  inplane  swirl  velocity,  9?,  averaged  over  incremen¬ 
tal  annular  regions  with  radii  (<^)  measured  from  the 
swirl  center.  The  upper  plot  in  the  right-most  column 
shows  an  enlargement  of  the  inplane  velocity  (ur+fi^r)) 
limited  to  a  region  twice  the  diameter  of  the  vortex 
core.  The  two  circles  identify  the  size  of  the  vortex 
core  (solid  curve)  and  the  vortex  void  (dashed  curve). 
Finally,  values  are  given  at  each  wake  age  for  the  void 
center,  the  vorticity  center  (or  centroid  of  vorticity), 
the  swirl  center,  the  vortex  convection  velocity  (ob¬ 
tained  from  the  path  of  the  vortex  voids),  and  the 
core  radius  along  with  the  maximum  swirl  velocity 
(both  of  which  are  fixed  by  the  peak  value  found  from 
the  curve  fit  io  ip  vs.  C).  These  values  are  also  given 
in  tables  3  to  5. 

Based  on  a  review  of  the  data  for  all  wake  ages 
(figs.  15  -  39),  the  highest  velocities  are  found  to  oc¬ 
cur  within  the  first  15  degrees  behind  the  rotor  blade. 
During  this  interval,  the  maximum  inplane  (swirl)  and 
out-of-plane  velocities  reached  44  and  12  percent  of 
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the  rotor  tip  velocity  (87.2  m/s),  respectively.  Al¬ 
though  the  degree  of  vortex  wander  was  somewhat 
erratic,  it  generally  increased  with  wake  age.  On  av¬ 
erage,  the  standard  deviation  (based  on  the  mean  lo¬ 
cation  of  the  void  center)  was  only  about  9  percent 
of  the  diameter  of  the  vortex  core  (the  worst  case 
being  18  percent).  The  calculated  locations  of  the 
vorticity,  swirl,  and  velocity-peak"*  centers  generally 
differed  from  the  void  center  by  about  11,  26,  and  11 
percent  of  the  vortex  diameter,  respectively  (fig.  40). 
There  was  no  apparent  correlation  between  the  void 
and  swirl  centers  from  image  to  image. 

The  out-of-plane  velocity  exhibited  the  most  un¬ 
usual  appearance.  Aside  from  the  boundary  layer 
wake  that  is  shed  from  the  trailing  edge  of  the  ro¬ 
tor  blade  and  is  most  distinct  when  ^  <  30®  (figs.  15 
-  17),  there  appear  to  be  two  concentrated  zones  of 
oppositely  directed  flow  near  the  center  of  the  vortex. 
Because  both  of  these  zones  are  generally  offset  from 
the  vortex  center,  the  bisecting  profiles  taken  along  z 
=  Constant  do  not  reveal  the  magnitude  of  the  “ex¬ 
cess”  and  “deficit”  velocities  that  contribute  to  the 
overall  structure  of  the  trailing  vortex. 

Even  though  the  convection  velocity  of  the  vor¬ 
tex  at  each  wake  age  has  been  subtracted,  the  rela¬ 
tive  “vertical”  component  of  velocity,  Wr^  exhibits  an 
asymmetric  profile.  Greater  symmetry  near  the  vor¬ 
tex  may  be  achieved  by  subtracting  the  local  down- 
wash,  which  can  only  be  approximated  from  a  time 
average  of  the  flow  field  at  corresponding  locations 
over  one  or  more  revolutions  (refs.  8,  10).  However, 
it  is  not  clear  that  this  effort  will  yield  any  better  un¬ 
derstanding  of  the  vortex  structure,  especially  because 
the  resulting  profiles  (including  the  peaJcs)  will  still  be 
asymmetric  (refs.  8,  13).  Given  the  inevitable  likeli¬ 
hood  of  asymmetry,  traditional  estimates  of  the  core 
size  and  peak  velocity  would  still  have  to  be  based  on 
an  average  between  the  two  halves  of  the  asymmetric 
profile,  and  this  is  effectively  what  is  already  provided 
in  the  more  general  annular  averages  that  produce  the 
(f  vs.  C  plot-  Nevertheless,  the  profile  plots  are  im¬ 
portant  because  they  retain  certain  nuances  that  can 
prove  helpful  in  explaining  other  unexpected  features 
in  the  flow. 

During  the  initial  development  of  the  trailing  vor¬ 
tex,  ^  <  30®,  the  vorticity  within  the  vortex  reached 
a  peak  value  very  near  the  center  of  the  vortex,  and 
rapidly  diminished  to  near  zero  at  a  distance  not 
far  beyond  the  core  diameter.  For  larger  wake  ages, 

^  >  45® ,  the  peak  values  of  vorticity  appear  to  reside 
in  an  annular  region  defined  by  the  size  of  the  void, 
yet  are  clearly  within  the  core  diameter.  Assuming 
these  measurements  are  valid,  this  variation  in  vor- 

^The  velocity-peak  center  is  based  on  the  locations  that  de¬ 
fine  dh  (see  fig.  41). 


ticity  would  imply  that  much  of  the  interior  region  of 
the  vortex  was  not  fully  in  solid-body  rotation.  Evi¬ 
dence  for  this  can  be  found  in  the  vertical  component 
of  velocity,  uJr,  which  shows  a  change  in  the  slope 
near  the  vortex  center  (see,  for  example,  the  profile 
at  =  255®  appearing  in  fig.  31).  This  peculiar 
departure  from  the  traditional  view  of  a  vortex  has 
appeared  in  other  rotor  vortex  data  (ref.  6,  where  the 
flow  was  seeded  with  incense  particles  and  measure¬ 
ments  made  with  a  laser  Doppler  velocimetry  (LDV) 
technique,  and  ref.  23). 

To  be  sure,  other  explanations  (seeding  being 
among  them)  for  the  unusual  velocity  and  vorticity 
profiles  over  the  interior  of  the  vortex  are  plausible. 
Although  centripetal  effects  may  have  been  a  factor, 
the  extent  to  which  the  velocity  within  the  core  was 
significantly  altered  is  unknown  because  the  size  of 
the  particles  was  not  measured.  Another  factor  may 
have  been  the  rival  between  legitimate,  but  barely  vis¬ 
ible,  small  particles  and  the  presence  of  sensor  noise, 
the  result  being  an  abundance  of  interpolated  vectors 
within  the  void  region.  Until  these  issues  can  be  more 
thoroughly  addressed,  it  is  uncertain  whether  the  data 
show  an  onset  of  vortex  instability;  reflect  a  difference 
in  particle  dynamics  because  of  a  nonuniform  distri¬ 
bution  of  particle  sizes  in  a  highly  accelerated  flow; 
or  are  simply  due  to  a  processing  limitation  arising 
from  a  marginal  scattering  of  light  from  very  small 
particles  that  resulted  in  incorrect  interpolations. 

The  location  with  the  greatest  number  of  inter¬ 
polations  always  occurred  within  the  void  region.  The 
number  of  interpolated  values  contributing  to  the  av¬ 
erage  at  these  “worst-case”  locations  ranged  from  16 
to  60  percent.  Although  far  from  ideal,  no  location 
consisted  solely  of  interpolated  values. 

To  address  a  possible  flaw  in  the  logic  used  to  av¬ 
erage  the  data  during  the  analysis  phase,  the  case  for 
V'  =  195®  was  reexamined.  First,  the  results  for  single 
images  (hence,  unaveraged)  were  individually  studied 
to  see  if  the  nonlinear  slope  in  the  core  velocity  might 
simply  have  been  an  artifact  of  the  averaging  pro¬ 
cess.  Although  the  unaveraged  velocity  profiles  were 
noticeably  more  “noisy,”  the  peculiar  change  in  slope 
near  the  center  of  the  vortex  was  found  to  be  present 
and  similar  to  the  conditionally  averaged  results.  As  a 
final  check,  a  different  criterion  was  considered  for  ob¬ 
taining  the  conditional  average.  In  this  case  the  more 
classical  approach  of  using  the  midpoint  between  the 
velocity  peaks  (highest  maximum  and  minimum  val¬ 
ues  occurring  during  horizontal  cuts  across  the  vortex) 
was  adopted.  Again,  the  variation  in  slope  across  the 
vortex  was  essentially  unchanged  from  that  observed 
when  the  conditional  averaging  criterion  was  based  on 
the  locations  of  the  voids.  In  fact,  the  same  conclu¬ 
sion  was  reached  after  considering  the  swirl  centers  as 
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the  averaging  criterion. 

Another  curious  feature  of  the  vortex  is  its 
slightly  elliptical  shape.  This  was  also  observed  in  the 
flow  pattern  constructed  from  LDV  measurements  in 
the  near  wake  of  a  hovering  rotor  (ref.  24).  A  minor 
amount  of  ellipticity  can  be  attributed  to  the  fact  that 
the  axis  of  the  vortex  is  not  perpendicular  to  the  plane 
of  the  light  sheet.  This  would  cause  the  vortex  to  ap¬ 
pear  stretched  in  the  radial  and  axial  directions  (rela¬ 
tive  to  the  rotor  coordinates).  The  amount  of  stretch¬ 
ing  in  both  directions  was  nearly  the  same  during  the 
early  wake  ages,  thereby  increasing  the  apparent  core 
size  by  about  0.06  percent  and  having  essentially  no 
effect  on  ellipticity.  As  the  vortex  aged,  the  apparent 
stretching  became  progressively  more  biased  in  the 
axial  direction  (because  the  radial  displacement  wets 
approaching  a  constant),  but  the  degree  of  ellipticity 
caused  by  the  angle  of  the  vortex  axis  was  still  less 
than  0.4  percent  after  one  revolution.  Judging  by  the 
amount  of  ellipticity  observed  in  the  data  (sometimes 
reaching  an  aspect  ratio  of  about  1.23),  and  by  the  al¬ 
ternating  orientation  of  the  major  axis,  the  cause  for 
this  elliptical  flow  pattern  cannot  be  due  to  a  chang¬ 
ing  angle  of  the  vortex  axis  as  the  vortex  ages  (unless 
small- amplitude  undulations  are  superimposed  on  the 
vortex  trajectory,  although  the  data  in  figure  40  do 
not  show  this). 

To  study  this  phenomenon  in  more  detail,  the  di¬ 
mensions  of  the  vortex  (major  and  minor  axes)  were 
equated  to  the  distances  between  the  velocity  peaks 
found  from  horizontal  and  vertical  cuts  across  the  cen¬ 
ter  of  the  vortex  (fig.  41).  The  peak-to-peak  distances 
associated  with  the  horizontal  and  vertical  cuts  were 
designated  by  dh  and  dy,  respectively,  and  these  are 
shown  for  wake  ages  ^  ~  0®  to  390®.  Using  red  and 
green  color  bars  to  indicate  the  size  of  d^  relative  to 
dvj  the  data  suggest  that  the  elliptical  pattern  is  ro¬ 
tating.  This  figure  also  shows  that  the  nonlinear  be¬ 
havior  of  the  velocity  profile  within  the  core  does  not 
depend  on  the  direction  of  the  cut  (horizontal  or  ver¬ 
tical)  or  the  direction  of  the  major  axis  of  the  ellipse 
(horizontal  or  vertical). 

According  to  the  location  of  the  void  center,  the 
radial  and  axial  coordinates  of  the  vortex  with  increas¬ 
ing  wake  age  are  known  (fig,  42).  Based  on  curve  fits 
to  this  data,  the  corresponding  convection  velocities 
were  calculated.  These  data  indicate  that  the  highest 
radial  velocity  occurred  immediately  after  the  trail¬ 
ing  vortex  left  the  trailing  edge  of  the  rotor.  This 
figure  shows  that  the  magnitude  of  the  radial  velocity 
rapidly  decreased  and  then  approached  an  asymptot¬ 
ically  diminishing  value  near  the  completion  of  one 
rotor  revolution.  The  behavior  of  the  axial  velocity 
component  was  quite  different.  The  magnitude  of  the 
axial  velocity  was  lowest  immediately  after  the  trail¬ 


ing  vortex  left  the  rotor,  but  then  increased  until  its 
progress  was  noticeably  impeded  after  the  passage  of 
the  following  rotor  blade.  The  velocities  may  have  a 
slightly  different  appearance  if  curve  fits  to  the  data 
are  separated  into  segments  before  and  after  passage 
of  the  second  blade  (?/j  =  180®).  The  wake  contraction 
can  also  be  viewed  in  the  v-z  plane  (fig.  43),  which 
shows  the  characteristic  “knee”  due  to  the  influence 
of  the  second  blade,  and  the  increasing  dominance  of 
the  axial  velocity  as  the  wake  aged. 

Based  on  the  radial  distribution  of  the  swirl  ve¬ 
locity  at  each  wake  age  (figs.  15  "  39),  the  peak  (or 
inflection)  in  the  curve  fit  was  used  to  define  the  core 
size  and  the  maximum  swirl  velocity.  The  initial  di¬ 
ameter  of  the  vortex  core  was  about  10  percent  of  the 
chord  of  the  rotor  blade  (c  =  104  mm)  and  more  than 
doubled  its  size  after  one  revolution  of  the  rotor  (fig. 
43).  The  maximum  swirl  velocity  was  highest  during 
the  initial  formation  of  the  trailing  vortex,  decreased 
rapidly  during  the  first  half  revolution,  and  settled 
on  a  nearly  constant  value  (on  average)  during  the 
IcLst  half  revolution  of  the  rotor.  It  is  difficult  to  say 
exactly  how  the  intensity  of  the  vortex  (which  is  re¬ 
lated  to  the  maximum  swirl  velocity)  changed  with 
wake  age  because  of  the  unavoidable  influence  of  the 
downwash  on  the  magnitude  of  the  swirl  velocity  (the 
vortex  is  embedded  in  a  nonuniform  velocity  field) . 

To  examine  more  closely  the  structure  of  the 
trailing  vortex,  contour  plots  of  the  vertical  compo¬ 
nent  of  velocity,  w^y  and  the  out-of-plane  vorticity, 
Qoy  are  shown  at  wake  ages  ^  =  15®  and  195®  along 
with  profiles  taken  from  several  closely  spaced  hori¬ 
zontal  cuts  across  the  vortex  (figs.  44  and  45).  The 
velocity  profiles  are  clearly  more  linear  within  the  core 
of  the  vortex  at  ^  =  15®  than  they  are  at  V’  =  195® 
(fig.  44).  The  corresponding  vorticity  profiles  exhibit 
a  single  peak  when  crossing  the  vortex  at  -0  15®, 

whereas  a  double  peak  is  evident  along  several  cross 
sections  of  the  flow  at  =  195®  (fig.  45). 

Another  point  of  interest  is  the  radial  distribu¬ 
tion  of  circulation  across  the  vortex.  Using  the  swirl 
velocity  to  calculate  the  circulation,  and  normalizing 
both  ip  and  T  by  values  existing  at  the  core  radius 
(Cc)j  results  for  ^  =  15®,  195®,  and  375®  are  shown  in 
figure  46.  As  the  distance  increased  from  the  center 
of  swirl  (("  =  0),  the  circulation  rapidly  increased  un¬ 
til  reaching  a  point  just  beyond  the  core  radius,  after 
which  the  rate  of  growth  decreased  significantly. 

Similar  to  the  selected  views  of  circulation,  the 
radial  distribution  of  vorticity  was  examined  during 
the  development  of  the  trailing  vortex.  However,  be¬ 
cause  the  vorticity  center  and  the  swirl  center  (as  well 
as  the  void  center)  data  were  generally  not  collocated, 
the  normalized  vorticity  profile  was  referenced  to  the 
vorticity  center.  The  peak  vorticity  occurred  at  the 
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center  («  =  0)  when  ^  =  15^  (fig.  47).  At  subse¬ 
quent  wake  ages  {ip  =  195®  and  375®),  the  vorticity 
in  the  center  decreased  and  the  peaks  progressively 
shifted  toward  the  core  diameter.  A  similar  change  in 
the  distribution  of  vorticity  was  reported  in  an  earlier 
fixed-wing  experimental  study  using  a  vorticity  meter 
(ref.  25).  As  discussed  earlier,  it  is  not  known  to  what 
extent  this  characteristic  might  be  due  to  potentially 
inaccurate  data  within  the  void  region.  Nevertheless, 
the  fact  that  this  radial  shift  in  the  vorticity  peak  has 
been  observed  while  using  other  measurement  tech¬ 
niques  (vorticity  meter  and  LDV)  suggests  that  the 
PIV  data  are  authentic  and  that  the  vortex  was  in¬ 
deed  becoming  unstable. 

The  self-similar  nature  of  the  swirl  velocity  and 
the  vorticity  for  the  entire  set  of  wake  ages  (^  =  0® 
to  390®)  are  shown  in  figure  48.  The  swirl  velocity 
was  made  dimensionless  by  its  value  at  the  core  ra¬ 
dius  {ip/(pc)i  and  the  radius  (measured  from  the  swirl 
center)  was  made  dimensionless  by  the  core  radius 
(C/Cc)-  Some  data  points  for  the  dimensionless  swirl, 
Ipy  exceed  1.0  because  the  maximum  swirl  velocity,  (pcy 
was  based  on  curve  fits  to  (p  vs,  (  at  each  wake  age 
(figs.  15  -  39).  The  dimensionless  form  of  a  popular 
vortex  model  (ref.  21) 


has  also  been  included  in  the  plot.  The  best  fit  to  the 
data  occurred  when  n  =  3  and  T^/Tc  =  1.26.  Ac¬ 
cording  to  this  model,  the  circulation  around  a  con¬ 
tour  placed  at  the  core  boundary,  Pc,  was  about  79 
percent  of  the  total  far-field  (or  large-radius)  value, 
r*.  Note  that  at  the  core  radius  ^  =  1,  so  that 

equation  (3)  requires  that  F^/Fc  =  2^/”  =  1.26  when 
n  =  3. 

The  self-similar  distribution  of  accumulated  vor¬ 
ticity  can  be  found  from  the  product  of  the  local  vor¬ 
ticity  (w^)  and  the  incremental  area  over  which  it  acts 
{ArAz) ,  and  then  summed  over  regions  defined  by  the 
radius  from  the  vorticity  center  (/c).  In  principle,  the 
result  should  be  equivalent  to  the  variation  of  circu¬ 
lation  with  radius,  provided  the  swirl  and  vorticity 
centers  are  coincident.  Within  the  core  (ac  <  Cc)j 
the  distribution  of  vorticity  is  quite  similar  and  is 
well  represented  by  the  model,  IpC.  Beyond  the  core 
radius  the  accumulated  vorticity  asymptotically  ap¬ 
proaches  a  constant  value  and  on  average  agrees  with 
timt  predictedjby  the  vortex  model,  which  states  that 
pC  ->  1,26  as  C  becomes  large. 

The  core-radius  circulation,  Fc,  was  observed  to 
increase  slightly  with  wake  age  (fig.  49).  Because 
some  vorticity  was  found  to  exist  beyond  the  core  ra¬ 
dius  at  all  wake  ages,  and  because  the  core  radius  gen¬ 


erally  increased  with  wake  age  (fig.  49),  it  is  probable 
that  additional  vorticity  was  gradually  being  gathered 
into  the  core.  The  total  far-field  vorticity  appears  to 
be  fairly  constant,  F*  =  1.42  m^/s.  The  growth  in 
the  size  of  the  vortex  core  can  be  approximated  by 
the  following  model  (ref.  26,  and  further  developed  in 
ref.  24): 


Cc  =  (Cco + (4) 

where  J  =  16  provides  the  best  fit  to  the  present  data 
(for  laminar  flow,  J  =  1).  Data  included  from  another 
experiment  (ref.  24,  based  on  LDV  measurements) 
was  reported  to  be  best  fit  by  the  model  when  S  = 
8.  It  should  be  noted  that  S  strongly  depends  on 
the  initial  core  size,  (^cO)  and  that  lower  values  of  (cO 
require  lower  values  of  S,  The  results  presented  in 
reference  24,  which  are  for  a  single-bladed  rotor,  show 
a  similar  range  of  core  sizes  and  growth  rate. 

Based  on  an  isolated  vortex  model  (refs.  22,  27) 
that  assumes  a  turbulence  mixing  length  that  is  pro¬ 
portional  to  the  vortex  radius,  the  ratio  of  the  core¬ 
radius  to  large-radius  circulation  should  be  fixed  by 
the  vortex  Reynolds  number,  F*/z/.  However,  the 
data  for  this  test  (F*/i/  ^  1.2  x  10^)  fall  above  the 
curve  and  are,  on  average,  aligned  with  a  laminar 
value  for  Fc/F*  (fig.  49).  This  would  imply  that  S  1 
rather  than  16  as  required  for  a  good  fit  to  the  Cc  ^5. 
7p  data,  and  would  suggest  that  either  the  models  are 
imprecise  or  that  the  core  radius  is  not  as  large  as 
measured. 

Just  as  the  core-radius  circulation  was  less  than 
the  large-radius  circulation  during  the  first  revolution 
of  the  rotor,  the  large-radius  circulation  is  expected 
to  be  less  than  the  maximum  bound  circulation  on 
the  blade.  Although  the  blade  circulation  was  not 
measured  in  this  test,  a  value  can  be  approximated 
from  blade-element  momentum  theory  (ref.  28).  Fora 
blade  with  constant  chord,  the  total  thrust  coefficient 
can  be  stated  in  terms  of  the  dimensionless  inflow,  A, 
as 

Ct  =  [  4X^rdr  (5) 

Jo 

and  the  sectional  value  in  terms  of  the  local  circula¬ 
tion,  F^,,  on  the  blade  as 


=  (5) 

From  the  above  two  equations  it  follows  that  the  in¬ 
flow  velocity  and  the  blade  circulation  (both  are  func¬ 
tions  of  r)  are  related  by 


g-Fj, 

QRc 


=  4A2 


(7) 
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The  local  inflow  velocity  for  hover  can  be  stated  as 
(ref.  28) 


(*) 

where  6  is  the  local  blade  angle.  Because  the  climb 
rate  in  this  test  was  insignificant,  hover  conditions  are 
assumed  to  apply.  The  change  in  blade  angle  along 
the  span  for  the  rotor  used  in  this  test  is  shown  in 
figure  50.  If  the  blade  angle,  ^(r),  is  known,  and  if 
it  is  assumed  that  Cia  0.1,  the  local  inflow  velocity 
predicted  by  equation  (8)  can  be  determined  (fig.  50). 
However,  in  order  to  more  easily  perform  the  integra¬ 
tion  given  in  equation  (5),  simple  analytical  expres¬ 
sions  for  the  inflow  in  the  form  of  A  =  Aor”  have  been 
proposed  (refs.  22,  28).  The  results  for  n  =  0  (con¬ 
stant  inflow),  n  0.5  (linear  circulation),  and  n  =  1 
(linear  inflow)  are  shown  in  figure  50.  Although  this 
power  function  has  served  well  to  illustrate  various 
performance  trends,  the  inflow  used  in  the  present 
analysis  was 


A  =  Aoysin^r  (9) 

From  equations  (5),  (7),  and  (9)  the  maximum  bound 
circulation  on  the  rotor  blade  is 


max 


(10) 


Recalling  that  Ct  =  0.00657  and  a  =  0.0691,  then 
equation  (10)  states  that  (rf,)mar  =  2.13  m^/s  for 
this  test.  To  summarize: 


^  =  0.79  (11) 

=  0.67  (12) 

implying  that  the  core-radius  circulation  is  only  about 
53  percent  of  the  maximum  bound  circulation  on  the 
rotor  blade. 


3.  The  boundary  layer  wake  that  was  shed  along 
the  span  of  the  rotor  blade  disappeared  after  = 
30°;  however,  two  concentrated  zones  of  oppositely 
directed  axial  flow  were  observed  near  the  center  of 
the  vortex  throughout  the  first  rotor  revolution. 

4.  The  void,  vorticity  (or  centroid  of  vorticity), 
and  the  swirl  centers  occurred  at  distinct  locations, 
but  were  never  separated  by  more  than  26  percent  of 
the  vortex  diameter. 

5.  The  axially  directed  velocity  profile,  relative  to 
the  convection  velocity  of  the  void  center,  was  asym¬ 
metric  because  of  the  influence  of  the  rotor  down  wash. 
The  vortex  also  appeared  to  have  an  elliptical  shape 
that  rotated  as  the  wake  aged. 

6.  The  peak  vorticity  occurred  near  the  center  of 
the  vortex  during  the  early  development  of  the  trailing 
vortex,  and  then  progressively  shifted  toward  the  core 
diameter.  It  is  presently  unknown  to  what  extent  this 
result  is  due  to  weakly  defined  characteristics  within 
the  void  region. 

7.  The  magnitude  of  the  radial  distribution  of  the 
circulation  rapidly  increased  until  reaching  a  point 
just  beyond  the  core  radius,  after  which  the  rate  of 
growth  decreased  significantly.  The  core-radius  circu¬ 
lation  increased  slightly  with  wake  age,  but  the  large- 
radius  circulation  appeared  to  remain  relatively  con¬ 
stant. 

8.  The  radial  distributions  of  swirl  velocity  and 
vorticity  exhibit  self-similar  behaviors,  especially 
within  the  core.  According  to  vortex  models  that  ap¬ 
proximate  the  measured  data,  the  core-radius  circu¬ 
lation  was  about  79  percent  of  the  large-radius  cir¬ 
culation,  and  the  large-radius  circulation  was  about 
67  percent  of  the  maximum  bound  circulation  that 
was  calculated  for  the  rotor  blade.  On  average,  about 
53  percent  of  the  maximum  bound  circulation  resides 
within  the  vortex  core  during  the  first  revolution  of 
the  rotor. 
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Table  3:  Locations  related  to  vortex 


Wake  age 
(deg) 

Vorticity  center 
r  (mm)  z  (mm) 

Swirl  center 
r  (mm)  z  (mm) 

Void  center 
r  (mm)  z  (mm) 

Core  radius 
(mm) 

Void  radius 
(mm) 

0 

945.78 

2.89 

946.11 

1.93 

946.00 

3.00 

5.00 

0.71 

15 

934.78 

-1.01 

934.25 

-1.31 

934.70 

-0.80 

5.25 

2.21 

30 

922.28 

-7.88 

922.99 

-8.66 

922.10 

-8.80 

5.00 

2.85 

45 

909.53 

-16.08 

909.51 

-14.88 

910.00 

-16.00 

6.75 

3.42 

60 

899.23 

-26.00 

899.35 

-25.80 

899.30 

-26.70 

6.75 

3.74 

75 

888.40 

-35.37 

888.98 

-34.31 

888.40 

-35.20 

6.75 

4.42 

90 

877.85 

-44.60 

878.04 

-45.51 

878.20 

-45.20 

7.50 

5.04 

105 

868.99 

-54.56 

867.06 

-54.17 

869.00 

-55.20 

7.75 

5.44 

120 

858.57 

-64.35 

858.66 

-63.02 

858.50 

-64.50 

8.50 

5.55 

135 

850.34 

-76.28 

851.04 

-73.69 

850.00 

-75.70 

8.50 

5.84 

150 

842.34 

-83.49 

839.95 

-84.16 

842.20 

-85.40 

10.50 

5.94 

180 

829.99 

-112.12 

828.45 

-112.83 

830.60 

-111.80 

8.75 

6.71 

195 

820.94 

-130.12 

820.29 

-130.54 

823.60 

-131.10 

8.75 

7.08 

210 

820.87 

-148.14 

819.50 

-146.08 

821.10 

-149.40 

9.51 

6.80 

225 

814.53 

-165.86 

812.56 

-166.17 

816.60 

-166.00 

8.75 

6.82 

240 

810.73 

-185.65 

808.70 

-183.76 

812.70 

-185.30 

9.99 

6.62 

255 

809.76 

-200.91 

808.47 

-201.46 

810.50 

-202.00 

10.24 

7.49 

270 

810.10 

-218.07 

807.03 

-217.92 

810.70 

-217.60 

9.75 

6.51 

285 

804.78 

-231.66 

801.58 

-231.76 

804.50 

-232.60 

10.00 

6.28 

300 

802.36 

-246.28 

800.52 

-245.89 

803.10 

-247.00 

9.75 

6.57 

315 

802.14 

-262.10 

800.28 

-261.70 

802.40 

-263.10 

9.75 

6.72 

328 

801.18 

-277.29 

803.30 

-276.32 

801.10 

-277.00 

10.00 

7.04 

344 

794.13 

-293.91 

791.73 

-295.04 

794.50 

-293.60 

10.50 

6.76 

375 

795.78 

-312.19 

796.70 

-309.02 

795.50 

-311.90 

11.00 

7.11 

390 

798.17 

-328.52 

797.79 

-325.25 

798.90 

-329.10 

11.25 

7.14 
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Table  4:  Vortex  convection  and  shape 


Wake  age 
(deg) 

Vortex  convection 
Uc  (m/s)  Wc  (m/s) 

Max  swirl  velocity 
(m/s) 

Velocity  peak  separation 
dh  (mm)  dy  (mm) 

0 

-3.43 

-3.03 

30.77 

10.77 

9.07 

15 

-3.43 

-3.04 

35.36 

9.58 

11.15 

30 

-3.41 

-3.07 

32.24 

10.93 

11.02 

45 

-3.37 

-3.14 

24.88 

13.64 

12.88 

60 

-3.24 

-3.31 

25.86 

13.51 

14.55 

75 

-3.23 

-3.39 

27.39 

13.68 

12.75 

90 

-3.12 

-3.57 

24.59 

13.65 

14.58 

105 

-2.99 

-3.79 

23.26 

13.67 

14.56 

120 

-2.84 

-4.04 

20.29 

13.49 

18.30 

135 

-2.66 

-4.30 

16.62 

14.91 

18.36 

150 

-2.48 

-4.57 

17,07 

16.26 

18.19 

180 

-2.09 

-5.07 

20.78 

17.61 

14.46 

195 

-1.89 

-5.28 

17.50 

19.03 

18.14 

210 

-1.70 

-5.43 

17.94 

18.04 

17.81 

225 

-1.52 

-5.53 

17.52 

17.99 

15.68 

240 

-1.35 

-5.58 

15.32 

19.37 

17.65 

255 

-1.20 

-5.58 

17.82 

19.38 

19.71 

270 

-1.06 

-5.53 

17.27 

19.02 

18.26 

285 

-0.95 

-5.45 

15.94 

18.99 

16.44 

300 

-0.85 

-5.35 

19.57 

18.97 

18.32 

315 

-0.77 

-5.25 

18.20 

19.04 

18.10 

328 

-0.72 

-5.17 

17.71 

19.07 

18.26 

344 

-0.68 

-5.08 

17.84 

19.02 

18.29 

375 

-0.64 

-4.99 

16.35 

18.96 

21.74 

390 

-0.64 

-4.99 

15.63 

19.11 

21.82 
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Table  5:  Vorticity  and  velocity  extrema 


Wake  age 
(deg) 

Vorticity  normal  to 
measurement  plane 
max(l/s)  min  (1/s) 

Velocity  within 
measurement  plane 
max  (m/s)  min  (m/s) 

Velocity  normal  to 
meeisurement  plane 
max  (m/s)  min  (m/s) 

0 

20463.3 

-2484.1 

36.32 

1.08 

4.79 

-10.78 

15 

22942.8 

-994.0 

38.76 

1.21 

2.62 

-9.40 

30 

21014.8 

-1290.1 

35.08 

1.54 

1.83 

-9.25 

45 

11374.9 

-1198.1 

26.59 

0.23 

2.33 

-7.61 

60 

11808.0 

-1055.6 

26.67 

0.02 

2.79 

-7.16 

75 

11929.1 

-1043.4 

29.03 

0.14 

2.71 

-6.72 

90 

10505.9 

-918.7 

26.52 

0.05 

4.57 

-6.61 

105 

9416.3 

-1028.0 

25.78 

0.05 

4.49 

—6.56 

120 

7657.5 

-1050.3 

23.93 

0.03 

4.69 

-6.55 

135 

5812.4 

-1172.7 

22.56 

0.20 

5.03 

-6.96 

150 

6990.6 

-1288.0 

25.69 

0.85 

4.55 

-9.30 

180 

8575.8 

-1143.9 

24.23 

1.37 

3.31 

-5.31 

195 

8755.1 

-1814.7 

23.55 

0.72 

3.52 

-6.09 

210 

8154.6 

-900.5 

23.14 

0.61 

3.48 

-5.76 

225 

7746.3 

-1147.3 

22.40 

0.44 

3.26 

-4.50 

240 

6852.9 

-1623.4 

19.61 

0.24 

3.46 

-4.41 

255 

7020.4 

-724.5 

20.62 

0.76 

4.19 

-5.36 

270 

7006.7 

-870.7 

22.13 

0.27 

2.77 

-4.82 

285 

6732.1 

-1318.6 

19.52 

0.06 

2.40 

-4.43 

300 

7062.0 

-706.5 

22.23 

1.37 

2.11 

-5.33 

315 

6382.8 

-773.3 

19.93 

1.20 

2.35 

-3.68 

328 

6532.7 

-788.4 

21.20 

0.90 

2.79 

-4.44 

344 

6841.8 

-859.0 

20.49 

0.38 

1.72 

-9.15 

375 

5851.2 

-492.8 

20.95 

0.73 

2.63 

-7.30 

390 

5977.2 

-734.8 

21.39 

0.05 

4.32 

-7.26 

Top  View 


Air  Exchanger  Section 


Figure  1:  The  7-  by  10-Foot  Subsonic  Wind  Tunnel  flow  circuit. 
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Top  View 


Flow 

Seeder 


Figure  2;  Instrument  setup  in  settling  chamber  of  7-  by  10-Foot  Subsonic  Wind  Tunnel. 
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Figure  3:  OLS  blade  profile. 


Figure  4:  OLS  model  rotor  head. 
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2Kx2K  120mJ 

CCD  Camera  Pulsed  Laser 

(85  mm  lens) 


Figure  6:  Particle  imaging  system  layout  (light  sheets  are  horizontal) . 
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Figure  7:  Perspective  effect  of  particle  displacement  with  stereoscopic  PIV. 


Figure  8:  A  single  interrogation  area  in  which  ten  particle  displacements  are  registered  on  a  22-  by  22-pixel  CCD 
grid. 


7^  X  magnification 
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Figure  9:  Particle  distributions  often  observed  around  apparent  vortex  voids  when  flow  is  not  uniformly  seeded 
(camera  view;  images  not  de warped). 
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Swirl  Velocity  Spline  Fit 


Adjust  for  Convection 


V  >  > 


Figure  11:  Sequence  for  extracting  vortex  characteristics. 


position 

vector 


a :  possible  swirl  center 
b :  neighboring  node 


Figure  12;  Angles  used  in  search  for  center  of  swirl. 
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Figure  13:  Wake  ages  at  which  PIV  data  were  recorded. 
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Figure  15:  Velocity  and  vorticity  characteristics  with  reference  to  the  center  of  swirl  ai  tp  ^  0^ , 
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Figure  17:  Velocity  and  vorticity  characteristics  with  reference  to  the  center  of  swirl  at  ^  30®. 
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Figure  18:  Velocity  and  vorticity  characteristics  with  reference  to  the  center  of  swirl  at  ^  =  45°. 
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Figure  19:  Velocity  and  vorticity  characteristics  with  reference  to  the  center  of  swirl  at  -0  =  60^. 
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Figure  23:  Velocity  and  vorticity  characteristics  with  reference  to  the  center  of  swirl  at  =  120®. 
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Figure  27:  Velocity  and  vorticity  characteristics  with  reference  to  the  center  of  swirl  3it  ip  =  195®. 
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Figure  29:  Velocity  and  vorticity  characteristics  with  reference  to  the  center  of  swirl  at  =  225^. 


Vortex  wander 


Central 

calculation 

grid 

Az=  2.0  mm 
Ar=  1.4  mm 


Case;  W11_255ES 
Measurement  Zone 


□ 


Collective=11°  Wake  age=  255° 

Void  center  r=  810.5  z=-202.0  mm 

Vorticity  center  r=  809.8  z=-200.9  mm 

Swirl  center  r=  808.5  z=-201.5  mm 

Vortex  conv  Uq=-1.20  Wc=-5.58  m/s 

Max  swirl  “  17.8  m/s 
Core  radius—  10.2  mm 
Void  radius=  7.5  mm 


-300 


30 


20 


10 


0  h 


-10 


.  vorticity  center 
^  apparent  void 


750 


800  850 

r,  mm 


900 


Detail  of  u  +  w^ 
within  ^  <  20.5  mm 
...  void  _  core 


10  20  30  40 

mm 


r,  mm 


Figure  31:  Velocity  and  vorticity  characteristics  with  reference  to  the  center  of  swirl  at  =  255®. 
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Figure  33:  Velocity  and  vorticity  characteristics  with  reference  to  the  center  of  swirl  at  =  285®. 
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Figure  35:  Velocity  and  vorticity  characteristics  with  reference  to  the  center  of  swirl  at  =  315®. 
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Figure  37:  Velocity  and  vorticity  characteristics  with  reference  to  the  center  of  swirl  at  =  344®. 
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Figure  38:  Velocity  and  vorticity  characteristics  with  reference  to  the  center  of  swirl  at  ^  =  375°. 
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Figure  39:  Velocity  and  vorticity  characteristics  with  reference  to  the  center  of  swirl  at  ^  =  390^. 
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Figure  40:  Void,  vorticity,  and  swirl  center  locations  from  ^  =  0^  to  390' 
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Figure  41:  Evidence  for  the  dynamic  elliptic  shape  of  the  vortex  during  ^  =  0^  to  390 


Figure  44:  Vertical  velocity  (rZir)  profiles  along  several  cross  sections  of  the  trailing  vortex  at  V'  =  15°  and  195°. 
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Figure  47:  Vorticity  (us),  calculated  along  radius  k  relative  to  the  vorticity  center  using  the  inplane  velocity  field  at 
=  195“,  and  375“. 
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Figure  48:  Nondimensional  swirl  velocity  and  vorticity  distributions  for  all  wake  ages  ip  =  0“  to  390” . 
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Figure  49:  Vortex  size  and  strength  during  wake  ages  ^  =  0^  to  390^. 
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